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With respect to usual thermal ferromagnetic transitions, the zero-temperature finite-disorder 
critical point of the Random-field Ising model (RFIM) has the peculiarity to involve some 'droplet' 
exponent 9 that enters the generalized hyperscaling relation 2 — a — vid — 9). In the present paper, 
to better understand the meaning of this droplet exponent 9 beyond its role in the thermodynamics, 
we discuss the statistics of low-energy excitations generated by an imposed single spin-flip with 
respect to the ground state, as well as the statistics of equilibrium avalanches i.e. the magnetization 
jumps that occur in the sequence of ground-states as a function of the external magnetic field. 
The droplet scaling theory predicts that the distribution dl/l 1+e of the linear-size I of low-energy 
excitations transforms into the distribution ds / s 1+e ^ di for the size s (number of spins) of excitations 
of fractal dimension df (s ~ l df ). In the non-mean-field region d < d c , droplets are compact df — d, 
whereas in the mean-field region d > d c , droplets have a fractal dimension df — 29 leading to the 
well-known mean-field result ds/s 3 ' 2 . Zero-field equilibrium avalanches are expected to display the 
same distribution ds/s 1+0 ^ d f . We also discuss the statistics of equilibrium avalanches integrated 
over the external field and finite-size behaviors. These expectations are checked numerically for the 
Dyson hierarchical version of the RFIM, where the droplet exponent 9(a) can be varied as a function 
of the effective long-range interaction J(r) ~ \/r d+a in d = 1. 



I. INTRODUCTION 

Since its introduction by Imry and Ma [l|, the random field Ising model where Ising spins Sj = ±1 interact 
ferromagnetically ( Jij > 0) and experience random fields hi of strength W 

E=~Y, JijSiSj -^hiSi (1) 

<i,j> i 

has remained one of the most studied disordered model over the years (see for instance the review Q and references 
therein). In addition to the usual critical exponents (a, ft, 7, v) which govern respectively the singularities of the 
free-energy f s i ng ~ |^| 2 ~ Q , magnetization m ~ \t\@ , susceptibility 7 ~ \t\ 7 , and correlation length £ ~ \t\~ v as in 
thermal ferromagnetic phase transitions, the zero-temperature phase transition that occurs at some finite disorder 
strength W c in the Random Field Ising Model (RFIM) has the peculiarity to involve a new exponent called 9 > 0, 
which enters the generalized hyperscaling relation (see for instance the recent discussion Q and references therein) 

2-a = v{d-0) (2) 

This specificity of the RFIM can be understood as follows. For usual thermal transitions, the free-energy singularity 
F(^ d ) associated to a correlated volume £ d is of order unity 



This corresponds to the density singularity 

T 
I 



F(Z d ) ~ T (3) 



fsing ~ja~\Tc- T\ d » ~ \T C - T\ 2 - a (4) 



and leads to the usual hyperscaling relation 2 — a = dv. For the RFIM however, the singular ground state energy 
E(£ d ) of a correlated volume £ d is not of order unity as in Eq. [31 but instead grows as some power of the size £ 

iW£ d ) ~ (5) 
so that the singular energy density is governed instead by 

e smg ~ |j ~ \W - W c \"^ ~\W- W c \ 2 - a (6) 
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leading to Eq. [2j As emphasized by Bray and Moore [H , the presence of 9 is directly related to the zero-temperature 
nature of the fixed point : for a thermal fixed point occurring at a finite T c , the fixed point corresponds to a fixed 
ratio Jl/T c so that the renormalized coupling Jl is invariant under a change of scale L; for the RFIM however, the 
fixed point correspond to a fixed ratio Jl/Wl where Wl represents the renormalized disorder strength. So even if 
this ratio is fixed, both the renormalized coupling Jl and the disorder strength Wl can actually both grow as L with 
the scale L. This explains how 9 enters the singular part of the energy in Eq. Another insight into the physical 
meaning of 9 comes from the following scaling picture : in a correlated volume of size £ d , the random fields correspond 
to an averaged external field of size 

1 t 

M^^Efti^r^ (7) 

which is expected to produce a magnetization of order 

m ef f{e)^xhefM d )^xC d ' 2 (8) 
and to lead to the following singularity in the energy density 

e sing ~ m eff (Z d )h eff (Z d ) ~ xC* *\W- W c \- 1+dv (9) 
If this analysis is correct, the exponent 9 of Eq. [5] should be directly related to the exponents 7 and v 



v 



(10) 



Since this analysis neglects possible correlations, Eq. [TU] is usually not considered as true, but it can be converted 
into the following bound 0, |j| 

e>- (11) 

v 

Taking into account the identity 2 — a = 2/3 + 7 and the generalized hyperscaling relation of Eq. [2l the inequality [10] 
can be rewritten as 

e>l- P - (12) 

2 v 

Recent works based on nonperturbative functional renormalization Q have concluded that the equality (|10[) does 
not hold in general. However numerically, it seems difficult to obtain a clear evidence, since the inequalities (|11[) or 
(|12p are actually satisfied as egalities within numerical errors for the short-range model in various dimensions : in 
dimension d = 3, the value of ^ turns out to be extremely small of the order S. ~ 0.012 Q, and the droplet exponent 
9 ~ 1.49 is not really optimized with respect to the non-optimized value d/2 = 3/2. In dimension d = 4, one has 
a 'reasonable' finite value ^ ~ 0.19 [1] and the droplet exponent 9 ~ 1.82 8i] remains close to d/2 — -. 

Besides this thermodynamic analysis, one expects that the exponent 9 has also the meaning of a 'droplet exponent' 
as a consequence of Eq. [5] In particular, it should govern the statistics of critical droplets, defined as low-energy 
excitations, as well as the statistics of 'equilibrium avalanches', i.e. the avalanches between ground states as a function 
of the external field H. The aim of the present paper is to study in details these statistical properties of droplets 
and avalanches. The various predictions are checked via numerics for the one-dimensional Dyson hierarchical version, 
where large systems can be studied, and where the value of the exponent 9 can be varied as a function of the exponent 
a of the long-range ferromagnetic interactions. Previous studies on these questions for the short-range model in 
dimension d = 3 can be found in 0, [TJj] for the statistics of low-energy excitations, and in [n], [TJJ for the statistics of 
equilibrium avalanches. 

The paper is organized as follows. After a reminder on the RFIM in the presence of long-range interactions in 
section [Til we describe the thermodynamical properties of the Dyson hierarchical model in section IIII1 In section IIV1 
we discuss the statistics of critical droplets, i.e. of low-energy excitations generated by an imposed single spin-flip 
with respect to the true ground-state. Section [V] is then devoted to the statistics of equilibrium avalanches, i.e. the 
magnetization jumps that occur in the ground state as the external field is varied. In section I VII we discuss the 
finite-size properties in the mean-field region d > d c where usual finite-size scaling is known to break down. Our 
conclusions are summarized in section I VIII In Appendix \K\ we describe how the sequence of ground states as a 
function of the external field can be constructed via an exact recursion for the Dyson hierarchical RFIM model. 
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II. REMINDER ON THE RFIM WITH LONG-RANGE INTERACTIONS 

For statistical-physics models in finite dimension d, a phase transition usually exists only above some lower critical 
dimension d > c?;. Thermodynamical critical exponents depend upon d in a region di < d < d c below the upper 
critical dimension d c , and then take their mean- field values for d > d c . For the short-range RFIM, one expects for 
instance di = 2 and d c = 6. Besides this short-range case, it is often convenient, both theoretically and numerically, 
to consider the case where the ferromagnetic coupling J^j in Eq. [T] is long-range and decays only as a power-law in 
the distance r — \i — j\ 

J(r) ~ ^ (13) 

where a > to have an extensive energy. Renormalization Group analysis of this long-range random-field model can 
be found in [H-[l3. 

A. Imry-Ma argument for the lower critical dimension di 

Let us consider the stability of the ferromagnetic ground state where all spins are (+), with respect to small random 
fields. If we flip a domain v ~ l d spins, the ferromagnetic cost scales as the double integral of Eq. [13] 

E DW {1) oc l d ~ a (14) 

whereas the random fields may correspond to an energy gain of order 

i d 

E RF (l)=Y / h l ocWli (15) 

z=l 

This argument yields that the lower critical dimension is 

d t = 2a (16) 

Indeed for d > di, E DW (I) > E RF (l) for sufficiently large ?, so the ferromagnetic ground state is stable with respect 
to small random fields, and one needs a finite critical disorder W c > to destroy it. On the contrary for d < di, 
E DW (I) < E RF (l) for sufficiently large I, so the ferromagnetic ground state is unstable with respect to small random 
fields, i.e. W c = 0. In the short-range case, the domain-wall cost scales as the surface E DW (I) oc l^ 1 yielding the 
usual value di = 2. 

B. Mean-field region d > d c 

The mean-field exponents for the thermodynamical observables at H = in the thermodynamic limit L — > +oo are 
given by the usual values 

e sing {W) ~\W- W c \ 2 ~ aMF with a MF = 
m{W) ~ \W - W c \ pMF with (3 MF = - 

X (W) ~\W- W C \-' MF with 7 mf = 1 (17) 
The Gaussian nature of the fixed point yields that the connected correlation (averaged over disorder) 



C connected (r) = < SoS r > — < So >< S r > (18) 

reads for a < 2 

C connected^) ~ / ^ C connected {<}) with Cconnected(l) ~ TT77 TTr i , i ,„ (19) 

J \W - Wc\ + \q\ a 
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where the term \q\ a represents the leading \ow-q singularity of the Fourier transform of the coupling J(r) ~ \ j r d + a 
(the usual short-range case corresponds the usual term q 2 recovered for a > 2). The term \ W — W c \ with power unity 
comes from the compatibility with the value Jmf = 1 f° r the susceptibility since 



X (W) = / d d rC connected (r) ~ / d d r J d d qe ttlr C connected (q) 

= d connected (q = 0)~\W-W c \- 1 (20) 
The large-r behavior of the connected correlation of Eq. [19] is thus the following power-law at criticality 



( jW=W c 
connected 



(r) „ J d d qe^L K r *-« = _L_ with w = 2 u (21) 



which leads to the following finite-size divergence of the susceptibility at W c 

Xl(W c ) = f d d rCZ=^ ed {r) = f d d rr^ d ~ L a (22) 

Off criticality, the exponential decay of the connected correlation defines the correlation length £ 

Cconnected(r) ~ / dqe^ ^ cx e~f with £ ~ |W - ^1"^'^ and vmf = - (23) 

The upper critical dimension d c is the dimension d where the generalized hyperscaling relation of Eq. [5] is satisfied 
with the mean-field exponents 

2 - a MF = VMF{d c - M f) (24) 
where the mean-field exponent 9mf is given by Eq. [10] 

8mf = ?ML = ^ (25 ) 
vmf vmf 



leading to 



d c = (26) 

vmf 



For a < 2 where vmf = 1/cr (Eq. |23|) . this yields 

d c (cr < 2) = 3(7 (27) 
whereas the usual short-range value d c — 6 is recovered for cr > 2 where vmf — 1/2 and 0mf = 2. 

III. DYSON HIERARCHICAL MODEL WITH RANDOM FIELDS 

To better understand the notion of phase transition in statistical physics, Dyson [17J has introduced long ago a 
hierarchical ferromagnetic spin model, which can be studied via exact renormalization for probability distributions. 
In this approach, the hierarchical couplings are chosen to mimic effective long-range power-law couplings in real space, 
so that phase transitions are possible already in one dimension. This type of hierarchical model has thus attracted 
a great interest in statistical physics, both among mathematicians (l8l - [2l| and among physicists [22Tj25| . In the field 
of quenched disordered models, Dy son hierarchical models have been introduced for spin systems with random fields 
[26[ or with random couplings |27l429l |. as well as for Anderson localization (30l - [37j . In the following, we consider the 
Dyson hierarchical random field model [26| . 



5 



A. Definition of the model 

The Hamiltonian for 2^ spins in an exterior magnetic field H reads 



2" 



n 2 n(H;{h 1 ,...,h 2 N}) =-J2( H + hi)Si (28) 

»=i 

- Ji [(Si + s 2 f + (S 3 + S 4 ) 2 + (S 5 + S 6 ) 2 + (S 7 + S 8 ) 2 + ...] 
- Ja [(Si + 5 2 + S 3 + Sif + (S 5 + S 6 + 5 7 + S 8 ) 2 + ...] 

- J 3 [(Si + 5 2 + S 3 + S 4 + S 5 + S 6 + S 7 + S 8 ) 2 + ...] - .... 

- Jn (Si + 5*2 + ... + S 2 n_i + S 2 n) 

The random fields hi represent quenched disordered variables drawn with some distribution, for instance the box 
distribution of width W 

, B „,„, 1 . ( w . , . w 



Pw(h) = W — < h < —j (29) 
The ferromagnetic couplings J„ are chosen to decay exponentially with the level n of the hierarchy 

To make the link with the physics of long-range one-dimensional models, it is convenient to consider that the sites i 
of the Dyson model are displayed on a one-dimensional lattice, with a lattice spacing unity. Then the site i = 1 is 
coupled via the coupling J„ to the sites 2™ _1 < i < 2™. At the scaling level, the hierarchical model is thus somewhat 
equivalent to the following power-law dependence in the real-space distance L n — 2™ 

J "={-hT" (31) 

One thus expects that the Dyson hierarchical model will have the same essential properties as the long-range model 
discussed in the previous section for the case d — 1. In particular, the model is well defined with an extensive energy 
for a > 0. The lower critical dimension di = 2a (Eq. I16j) coming from the Imry-Ma argument and the upper critical 
dimension d c = 3er (Eq. |2"T)) coming from the analysis of the Gaussian fixed point are the same : so in d = 1, the 
mean-field region and the non-mean-field region exist respectively in the following domains of the parameter a 

mean — field region d = 1 > d c for < a < — 

3 

non — mean — field region di < d = 1 < d c for - < a < ^ (32) 



B. Value of the droplet exponent 8 = a in the non-mean-field region 

Whereas for the short-range model, we are not aware of any conjecture concerning the value of 9 in the non-mean- 
field region di = 2 < d < d c = 6, Grinstein j 1 31 ] has conjectured that in the presence of long-range interactions, the 
exponent 8 always takes the following simple value 

6 = a (33) 

which is consistent with various limits (in particular d — > di = 2a and d — > d c = 3a) and various perturbative 
expansions jl3-15]. This conjecture was then shown to be wrong perturbatively at order 0(e 2 ) [l4[ for the long- 
ranged case discussed in section [Til but to be true in the Dyson hierarchical version [26[ . 

Eq[33]can be understood via the following scaling analysis. Let us introduce the exponent y governing the power-law 
decay of the magnetization per spin with the system size L exactly at the critical point W c 



m(W c ,L) cx L~ y 

L— y-\-oo 



(34) 



6 



The effective magnetic field per spin resulting from the random fields scales as 

^h iL * L-i (35) 

i=l 



h e / f 



The characteristic Zeeman energy of the L d spins of magnetization Ml = L d m(W c , L) is thus 

E z (W c ,L)~h e L ff M L oc Li-y (36) 

whereas the ferromagnetic energy associated to the effective coupling at scale L (see Eq. ITU]) 

rjf „ L- d - a (37) 

behaves as 

E ferro (W ei L)~Jl ff Ml oc L d -°- 2 y (38) 

L~ > +oo 

At the critical point, the two energies of Eq. [55] and Eq [33] should remain in competition at all scales, i.e. they 
should have the same scaling 

y = \-° (39) 

Then, the two energies scale as 

Ez(W c ,L) ~ E ferro (W c ,L) ~ L e with 6 = a (40) 

So Eq. [33] relies on the fact that the renormalized ferromagnetic coupling is directly given by the power-law 

defining the model (Eq. 1371) as a consequence of the exact hierarchical structure. As a final remark, let us mention 
that in the short-range case, the effective renormalized ferromagnetic coupling cannot be simply estimated, and 
this is why there is no simple conjecture for the value of 9. 



C. Finite-size properties of thermodynamical observables 

As a consequence of Eq. |33]that holds exactly for the Dyson hierarchical model, many finite-size behaviors exactly 
at criticality can be explicitly computed as a function of a, for all d > di, i.e. both in the non-mean field region 
d < d c and in the mean-field region d > d c (since in the mean-field region, one has also 9mf = & fEqs 1231 and 1251 )). 
In particular, the singular part of the energy density scales as (Eq. I40[) 

e smg (W c ,L)~ — [ =L a - d (41) 
the divergence of the susceptibility scales as (Eqs [3411351 and [39]) 

X (W c ,L)~?^^~L d />-*=L° (42) 

The exponent of the two-point correlation may also be obtained as 

SA7r~ m 2 (W c ,r) ~ -L = with f, = 2 - 2a (43) 

At low temperature T, the connected correlation involves another exponent 

< S S r > - < So >< S r > ~ (T) r - {d - 2+ ^ (44) 

As a consequence of the scaling of Eq. 001 only a rare fraction T/r 9 can contribute to the connected correlation 
function, so that one expects the following shift 

V = v + 8 = 2-<7 (45) 
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D. Numerical results on the magnetization 

As explained in Appendix the hierarchical structure of the Dyson model allows to write exact recursions to 
compute the ground states that occur as a function of the external field for each disordered sample. We have studied 
the following sizes 2 6 < L < 2 21 with a corresponding statistics of 4.10 7 > n s (L) > 45. 10 3 independent samples. 



1. Location of the critical point 




FIG. 1: Location of the critical point W c via the crossing of the rescaled magnetization curves L 1,/2_tT \m(W, L)\ as a function 
of the disorder strength W, for sizes 2 6 < L < 2 21 . (a) Case a = 0.4 (Non-mean-field region) : the curves L 01 \m(W, L)\ cross 
near W c c± 8.8 (b) Case a = 0.2 (Mean-field region): the curves L 3 \m(W, L)\ cross near W c m 23. 

Exactly at criticality, the magnetization |m(W c , L)\ is expected to decay as L~ y with y = 1/2 — a (see Eqs [34l and 
I39p . both in the mean- field region and in the non-mean- field region. With our numerical data, we indeed find that the 
curves L 1 / 2- CT m(W, L) for various sizes L cross more and more sharply as L grows : this crossing allows to locate the 
critical disorder strength W c , as shown on Fig[T]for the two cases a = 0.2 and a — 0.4. We have also data concerning 
the cases a = 0.1 and a = 0.3 (not shown). All numerical results presented below concern the critical point W c , i.e. 
with our data 



W c (a 


= 0.1) 


~ 50 


W c (a 


= 0.2) 


~ 23 


W c {a 


= 0.3) 


~ 13.5 


W c (a 


= 0.4) 


~ 8.8 



'Non-mean-field' region | < a < i with usual finite-size scaling 



In the 'non-mean-field' region, thermodynamical observables follow standard finite-size scaling forms in (W — 
W C )L 1 / 1 ' involving the correlation length exponent v, so that the known behaviors in L at criticality for the magne- 
tization, susceptibility and singular energy yield 

P 1 

— = V = c 

v y 2 



2 

V 

2- a 
v 



(47) 
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4 

1/2- a 




FIG. 2: Finite-size scaling of the magnetization at W c as a function of the external mag netic field H : l} 12 -" m{H, L) as a 
function of the scaling variable HL 1 ^ 2 (see Eq. I48[) for a = 0.4 (non-mean-field region). 

Since the correlation length v is not known exactly, the values of (a, j3, 7) are not known either, but the ratios of Eq. 
l47l are exactly known. 

Besides these good finite-size properties in (W — W C )L 1 / V at H — 0, one also expects good finite-size properties in 
the variable HL 1 / 2 (see Eq. 155)) exactly at W c The following finite-size scaling form 

m(W c ,H> 0) oc L a - 1/2 M(HL^ (48) 

is indeed well-satisfied as shown on Fig. [5] for the case a = 0.4. 

The finite-size properties in the mean-field region d > d c where usual finite-size scaling is known to break down will 
be discussed later in section IVII 

In conclusion of this section, the Dyson hierarchical RFIM is a convenient model where large systems can be studied 
numerically with a good statistics, and where the droplet exponent 9 = a can be varied by choosing the parameter a 
of the effective long-range interactions. Now that we have located the zero-temperature finite-disorder critical point 
for various cr, we may study the statistics of low-energy excitations and of equilibrium avalanches in the remaining 
sections. 



IV. STATISTICS OF LOW-ENERGY EXCITATIONS AT CRITIC ALITY 

In disordered systems, there can be states that have an energy very close to the ground state energy but which are 
very different from the ground state in configuration space. In the droplet theory, developed initially for the spin-glass 
phase [38j and then for the frozen phase of the directed polymer in a random medium [39j - |4lj |. the low-temperature 
physics is described in terms of rare regions with nearly degenerate excitations which appear with a probability that 
decays with a power-law of their size. In these models of spin-glasses or directed polymers, the droplet exponent 9 is 
a property of the low-temperature disorder-dominated phase T < T c . In the RFIM however, the originality is that 
the droplet exponent 9 is a property of the zero-temperature finite-disorder critical point. We have already recalled in 
the introduction how this exponent 9 enters in the critical properties of thermodynamical observables, in particular in 
the generalized hyperscaling relation of Eq. [5] In the present section, we discuss how this droplet exponent 9 governs 
the power-law distribution of low-energy excitations at W c . 

A. Low-energy excitations generated by an imposed spin-flip 

To generate low-energy excitations above the ground state in disordered systems, various procedures have been 
followed in the literature (see for instance the various methods concerning spin-glasses |43 - [4q ). In the following, 
since we wish to generate an elementary local excitation, we have chosen the 'single spin flip method', already used 
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in Ref. [Hj for the short-range 3D RFIM. The idea is the following : in each disordered sample, we first compute 
the true ground-state. Then we impose the flip of a given spin Sj with respect to its orientation in the true ground 
state, and we compute the new modified ground state when this constraint is taken into account. We measure the 
number s of spins that are different in the two ground states. This excitation has a finite-energy by construction : if 
only Si flips, the cost is simply AE = 2\h\° c \ in terms of the local field h\° c ; if the systems chooses to flip s spins, it 
is because the energy cost is lower. 

At criticality, the probability distribution (s) of the number s of spins of this low-energy excitation is expected 
to decay as a power-law in the thermodynamic limit L — > +00 

DY c (s) ~ — (49) 

where td is defined by this equation. For finite L, a finite cut-off 



s*(X) cx L p (50) 

L— ^-+00 



is expected to govern the far-exponential decay 



D?°(s) * e-^m (51) 

s— > +00 

Let us first discuss the relation between the exponent td of Eq. HHl and the droplet exponent 9. 



B. Relation with the statistics of the linear size I of droplets 



From the definition of the droplet exponent 8, a droplet of linear size I has an energy cost of order 

E dr °P let (l) cx l e u (52) 

where u is a positive random variable of order 0(1) distributed with some law p{u) . The random variable u is expected 
to have a zero weight p(u = 0) > at the origin. The probability to have a droplet of linear size I and of finite energy 
E droplet (l) < Eq then scales as Prob(u < ~ p(u = Taking into account the logarithmic measure dl/l in the 

size / to insure independent droplets, one obtains the following distribution for the linear size I 

dlP *ro P let {l) ^_M_ Eop{u ^ 0) (53) 

as in other models like spin-glasses [38| or directed polymers (39rl4l1 |. 

To obtain the statistics in the size s (number of spins) of droplets, one needs to know the fractal dimension df of 
droplets 

s - l d t (54) 

The probability distribution in / of Eq. [53] transforms into the following power- law distribution of Eq. [49] with 

, , 

td = 1 + -r (55) 
d f 

In the context of spin-glasses, a similar relation between the droplet exponent 9 and the avalanche exponent r 
(which coincides with the exponent of low-energy excitations of Eq. [53] as discussed in the next section) has been 
discussed in Ref. [57}, where it is assumed that the density of droplets is (dl/l 1+e ) x l d ~ d f which is different from 
Eq. [53] whenever the droplets are non-compact df ^ d. We believe that Eq. [53] is correct, and is consistent with our 
numerical results where we generate droplets by flipping an arbitrary spin , both in the case of compact or fractal 
droplets, as we now discuss. 



C. Statistics of the size s of low-energy excitations in the non-mean-field region d < d c 



In the non- mean- field region, one expects that droplets are compact, i.e. the fractal dimension of Eq. [M]is simply 

d f = d (56) 
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In D(s) 





FIG. 3: Distribution D£ c (s) of the size s (number of spins) of critical droplet for a — 0.4 (non-mean- field region) : (a) ]nDf"(s) 
as a function of Ins to measure the droplet exponent td in Eq. 1491 : the slope corresponds to the exponent td = 1 + = 1.4 
(b) Finite-size scaling In (L pTD D^ c (s)) as a function of In (^) with p = 2a = 0.8. 



so that Eq. l55l reads 

Q 

t d = 1 + - for d < d c (57) 
a 

By considering avalanches in the next section \V\ (see Eq. l91"j) . we expect that the cut-off exponent p is directly related 
to the droplet exponent 9 

p = 26 (58) 

Our numerical data for the Dyson hierarchical Dyson model of parameter a — 0.4 are in agreement with these 
prediction. On Fig. [3] (a), we show a log- log plot of the droplet distribution D^ c {s) for various L and we measure 
the slope td(<J — 0.4) = 1 + 6 = 1 + a = 1.4. On FigO (b), we show that a satisfactory data collapse can be obtained 
in terms of the rescaled variable with p(a = 0.4) = 29 = 2a = 0.8. 



D. Statistics of the size s of low-energy excitations in the mean- field region d > d c 

In the mean-held region d > d Cl one expects that 'loops' are not important, so that the probability distribution 
D(s) to return s spins to obtain the new ground state when one imposes the flip of one given spin Si should take 
the same form as the probability that a spin belongs to a cluster of size s for the percolation problem of the Bethe 
Lattice that has no loops [i?} 

D M f(s) = J— e s(w-w c f with tmf = I (59) 
This means that the density n(s) of clusters of size s per spin reads 

n(s) = = _J_ e -s(w-w c f m 

Let us now adapt the Coniglio scaling analysis [481 ] concerning the percolation transition for d > d c to our present case. 
In a volume l d , the number N[d of clusters has the following singularity in terms of the cut-off s max ~ (W — W c )~ 2 
of Eq. M 

N ld =l d J2 n(s) = l d [s max r™ F = l d (W - W c f T ^ (61) 
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FIG. 4: Distribution D™" (s) of the size s (number of spins) of critical droplet for a — 0.2 (mean-field region) : (a) In D™ a (s) 
as a function of Ins to measure the droplet exponent td in Eq. [49] : the slope corresponds to the mean-field exponent 

to = t d if = 3/2 (b) Finite-size scaling In (l pt " " D™'" (s)) as a function of In with p = 2a = 0.4. 



In particular, on a correlated volume £ rf = \W — W c \ " MF , the number of clusters has for singularity 

N id = £ d (W - W C ) 2TMF = £, d ~ dc with d c = = — (62) 

vmf vmf 

As stressed by Coniglio [48| , the meaning of the upper critical dimension d c is that it separates : 

(i) the mean-field region d > d c where the number N^d of clusters diverges with £ as £ d ~ d <= 

(ii) the non-mean-field region d < d c where a correlated volume contains a single relevant cluster. 

In our present long-range case where vmf = 1/(7 (Eq. [23]) we recover d c = 3a (Eq [27]). whereas the short-range 
case where ^/i? = 1/2 corresponds to d c = 6. 

In the mean-field region, the fact that a correlated volume contains a large number of clusters Ned 3> 1 is possible 
because each cluster has a fractal dimension d/ given by 

W ~ (W - VK C )- 2 - e d/ with d f = — = 26 MF (63) 

"MF 

in agreement with the relation of Eq. [55] 

r MF = l + — = - (64) 
In our present long-range case where vmf = l/* 7 (Eq. I23p . the fractal dimension of droplets is thus 

d f = = 2cr (65) 

vmf 

whereas in the short-range case where vmf = 1/2, the fractal dimension is df = 4. 

In this mean-field region where droplets are non-compact, the only constraint on the maximal linear size of droplets 
is the system-size L itself : l m ax ~ L, so we may expect that the maximal size in s of avalanches in a finite system is 

s*(£) - L d f = L 2a (66) 

Our conclusion is that in the mean- field region, the exponent p of the finite-size cut-off of Eq. [50] should be 

p = d f = 2a (67) 

Our data for a — 0.1, a — 0.2 and a = 0.3 are indeed compatible with the values tmf = 3/2 and p = 2a, as shown 
on Fig. g]for a = 0.2. 
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V. STATISTICS OF EQUILIBRIUM AVALANCHES 



Power-law avalanches occur in many domains of physics. In the field of disordered systems, non- equilibrium 
avalanches have been much studied, in particular in the context of driven elastic manifolds in random media (see the 
review [Jj| and references therein, as well as the more recent works [13, HH) and in the Random Field Ising model (see 
[52l - [54j and references therein). To better understand the properties of these non-equilibrium avalanches, it seems 
useful to make the comparison with the equilibrium avalanches that occur in the same systems. The statistics of 
equilibrium avalanches has been studied for elastic manifolds in random media [56], for spin-glasses [57| and for the 
RFIM [ll|, Oj)] , where some universality has been found between equilibrium and non-equilibrium avalanches [lH . In 
this section, we discuss the statistics of equilibrium avalanches in the RFIM. 



A. Observables concerning equilibrium avalanches 



Let us introduce the average number of avalanches of size s that occur at H in a system of size L 



N L {s,H) = Y,5{H-h }llp ( l )) S(s- mi+1 m ) (68) 



where the hfn p (i) and the rrii are the fields and the magnetization occurring in the sequence of ground-states of a 
sample as a function of the external field (see more details in Appendix [a"| . The total number of avalanches for 
— 00 < H < +00 is then 



/+00 p+ao 
dH / dsN L (s,H) 
-00 Jo 



Nl ot = dH dsN L (s,H) (69) 



Since during the history, each spin of the L d spins has to flip exactly once, one has the exact sum rule 

/+OO /- + 00 

dH I dssN L (s,H) = L d 
-00 Jo 

In terms of these avalanches, the difference between the magnetization at H and at H = can be written as 



so that the susceptibility reads 



dm 1 dM 1 



+ OO 



(70) 



rH p-\-oo 

M(H) - M(0) = 2 / dh dssN L {s,h) (71) 
Jo Jo 



^ H ^dH-17dH-I?J Q ^ 



B. Probability distribution of zero-field avalanches at W c 

To measure the probability distribution of the size s (number of spins) of avalanches exactly at H = (see the 
definition of Eq. EH]) 



>(i?=o)^ _ N L (s,H = 0) 

J +o ° dsN L (s,H = 0) 

we have considered, in each disordered sample, the avalanche occuring at the smallest \hfu p (i)\. 

h?™ = min\h fli M (74) 

The probability distribution R(hi oc ) of the local fields hi oc {i) > of spins Si in the ground state is expected to 
have a finite weight R(hi oc = 0) > at the origin hi oc = 0. Then the closest avalanche from H — occurs at a field 
of order 

- IS (75) 
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FIG. 5: Probability distribution P ( L H - U> (s) of the size s of avalanches at W c and H = in the non-mean-field region for a = 0.4 
and the sizes (2 6 < L < 2 21 ) : (a) \nP i j H ^°\s) as a function of Ins : the slope corresponds to r = 1 + a = 1.4 (b) Finite-size 
scaling analysis : In ( L pT Pf I ^°\s) J as a function of In (jj) with r = 1.4 and p = 2a = 0.8. 



| 1 1 1 1 1 1 1 10 




1 2 3 4 5 6 7 

In s 



-6 -5 -4-3-2-10 1 2 3 

ln(s/L P ) 



FIG. 6: Probability distribution P [ L ' (s) of the size s of avalanches at W c and H = in the mean-field region for a = 0.2 
and the sizes (2 6 < L < 2 21 ) : (a) In P^ H ^°\s) as a function of Ins : the slope corresponds to the exponent r = ta/f = f (b) 
Finite-size scaling of the probability distribution lnP£ H ~°^(s) of avalanches : In ^L pT p£ H ~ ' (s)^ as a function of In (-^-) with 
r = | and p = 2a = 0.4 . 

(If one draws L d variables from this distribution, the minimal local field will scale as h™™ oc 1/L d from the estimate 

! ° C dhR(h) = R(0)h™ n (76) 

We have checked the validity of Eq. [75] for the Dyson hierarchical RFIM, both in the non-mean-field and in the 
mean- field regions. 

When the external field H reaches this h^™ oc 1/L d , at least one spin becomes unstable, and it induces an avalanche 
of size s with some probability Pl(s). From this argument, it seems natural to expect that this Pl(s) exactly coincides 
with the droplet distribution -Dl(s) discussed in the previous section HVl (the only difference is that for the droplet 
distribution, we have chosen to force the flipping of an arbitrary spin, whereas here the field forces the flipping of the 




14 



spin having the lowest local held) 

P { L H=0) (s) =D t (s) = -U^^tc> with t = t d (77) 

with the same finite-size cut-off s*(L) ~ L p . As explained in section llVl the difference between the mean-field and 
the non-mean-field regions is that 

Q 

r(d< d c ) =T D {d<d c ) = l + - 

d 

r(d>d c ) =t% f = \ (78) 

Our numerical data for the Dyson hierarchical RFIM of parameter a = 0.4 (non-mean-field region) and a = 0.2 
(mean- field region) are in agreement with this picture, as shown on Figures [S] and [5] respectively. 

C. Integrated probability distribution of avalanches for —oo<H< +oo at W c 




FIG. 7: Integrated distribution P}J lt (s) of the size s of avalanches for — oo < H < +oo at W c in the non-mean-field region for 
a = 0.4 (sizes 2 6 < L < 2 21 ) (a) \nPl nt (s) as a function of Ins : the slope corresponds to the exponent Tint = f + a = 1.9 (b) 
Finite-size scaling analysis : In (L pTi " f Pj, (s)J as a function of In [jj) with T; nt = 1.9 and p — 2a = 0.8. 



We have also measured numerically the so-called 'integrated' distribution of the size s of avalanches (where 'inte- 
grated' means 'integrated over the external field — oo < H < +oo'), 

f + °° dHN L (s,H) 

PlKs) = J -°° Nr yi (79) 

where the notations Nl(s,H) and 7V£ ot have been introduced in Eqs IBS] and [S3 

Lets us now focus on the non-mean-field region d < d c , where one expects that the statistics of avalanches of Eq. 
[68] will follow some finite-size scaling form 

dHdsN L (s, H) ~ dHL a ^<P (JL. Hs^ (80) 

The exponents r and p have been already introduced as the exponents that characterize the statistics of zero-field 
avalanches of Eq. [73] 

(H=0) = N L (s,H = 0) = N L ( S ,H = 0) J_ / _s_ x 

[ >~ f Q + °°d S N L (s,H = 0) f Q + °°dsN L (s,H = 0) f*\lJ»") ^ 
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The exponent a governs the total number Nj"* of avalanches (Eq. l69l) occuring in a sample of size L 

dH J dsN L {s,H) ~ L a (82) 

The exponent ip describes the finite-size scaling in the external field H . In this non-mean-field region where avalanches 
are compact s ~ l d , we expect that the appropriate scaling variable is Hl d l 2 — Hs 1 ^ 2 so that ip takes the simple value 

tp = - for d < d c (83) 

The integrated distribution of avalanches introduced in Eq. [721 can be now written in terms of the finite-size scaling 
form of Eq. [SO] 

r + °° dHN L (s,H) i r + °° / s ,\ i r +ca / s \ 

T'W = ^ = F /_ (£, - ^ d» (£,*) (84, 

In the thermodynamic limit L — > +oo, the integrated distribution is also a power-law 

1 



Pri +00 (s)~— (85) 



r+oo 



s 

where the exponent Tint is shifted from the exponent r = 1 + g of the zero-field avalanches (Eq. 177)) by the factor 

^nonMF = l/2 (Eq.ES) 

3 ^ 

Tint = T + V = 2 + 3 ( 86 ) 

Our numerical data for cr = 0.4 are in agreement with this value Ti nt = I + a = 1.9 as shown on Fig. [JJ (a). 

Let us now write the consistency equations that fixes the values of p and a. Eq. [72] concerning the susceptibility 
yields 

XL ( H = ^ = T d ~J ds sN L (s,H = 0)=L a - d daa 1 -**^— ,0) ~i - d (i") (87) 

From the divergence of the susceptibility at criticality %l ~ i" = L 9 (Eq. [T0| . we obtain the relation 

6 = a-d + p{2-T) (88) 
The other consistency relation comes from the sum rule of Eq. [TD] that reads 

/+oo r+oo r+oo r+oo 

dH ds sN L (s,H) = L a dH ds s 1 " T $ ( — ,Hs 4 ') 

-oo Jo J-oo JO \LP J 

dH^h)~L a {LPf- T ^ (89) 

i.e. one obtains the following relation between exponents 

d = a + p{2 - t - ip) (90) 
The difference with the previous relation of Eq. [88] yields 

P = 7 = 20 (91) 

ip 

in agreement with the value measured for droplets (see Fig. [3](b)), for zero-field avalanches (see Fig. E](b)), and for 
integrated avalanches (see Fig. [3(b)). 

The exponent a governing the total number of avalanches (Eq. [82]) then reads from Eq. [88j Eq. [91] and r = 1 + ^ 

a = 6 + d- p(2-t) = d-8 + 2— (92) 

d 

Besides our numerical checks concerning the Dyson hierarchical RFIM, we may also compare with the numerical 
results of Ref. [12] concerning the statistics of equilibrium avalanches in the short-range RFIM in d = 3 : using the 
value of droplet exponent 9 ~ 1.49 0], we expect that zero-field avalanches correspond to the exponent r = 1 + g — 1.5 
(unfortunately very close to the mean- field value !), and that the integrated avalanches correspond to the exponent 
T int — 3/2 + 2 — 2-, in agreement with the value measured in Ref. 

In the non-mean-field region d > d c where usual finite-size scaling forms breaks down, the analysis of finite-size 
properties requires a more subtle analysis, as discussed in the next section. 
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VI. FINITE-SIZE PROPERTIES IN THE 'MEAN-FIELD' REGION d > d c 

As is well known, the usual finite-size scaling properties are not valid in the mean-held region d > d c (see for instance 
[58U62j and references therein). In particular, the equalities of Eq. H71 concerning thermodynamic observables are not 
valid anymore, because the correlation length £ ~ (W — W C )~ I/ is not the only important divergent length scale in 
the system. These properties can be more clearly understood by the scaling theory developed by Coniglio [48| on the 
example of the percolation transition for d > d c . In section HV Dl concerning the statistics of low-energy excitations 
in the mean-held region d > d c , we have already started to explain how Coniglio's approach could be adapted to the 
RFIM. In the present section, we continue this analysis and derive the consequences for the finite-size properties of 
various observables. 



A. Finite-size properties at H = as a function of (W — W c ) 

In section HV Dl concerning the statistics of low-energy excitations in the mean- field region d > d c , we have recalled 
how the irrelevance of loops for d > d c leads to Eq. [SD] for the density of droplets, and to the fractal dimension 
df = 2/vmf = 29m f (Eq. H3]). We have already described how the number of clusters N^d in a correlated volume £ d 
grows as £ d_rf <= (Eq. |62|) . This means that the correlation length £ is not the only important length scale (in contrast 
to the region d < d c ). In particular, it is clear that another important scale is the smaller length £i dehned by the 
requirement that the number N^a of clusters in a volume £f is one 

% = £i E n ( fi ) = &~ da = 1 y iclds & = e* = (w- w c )- vm *% (93) 

s 

For instance in the short-range case where vmf = 1/2 and d c = 6, one has £ ~ (W — Wc)^ 1 ^ 2 and £i ~ (W — W c )~ 3 / d . 
In the one-dimensional (d = 1) long-range case where Vmf = 1/<t and d c — 3cr, one has £ ~ (W - Wc)- 1 '* and 
£i ~ (W-W c )~ 3 . 

In a finite system of linear size L, one may thus expect three regimes : 

(i) for L < £i : there exists N C (L) = 1 cluster of singular energy E C (L) = L e 

(ii) for £i < L < £ : there exists N C (L) = ^ clusters of singular energy E C (L) = L 8 

d 

(hi) for £ < L : there exists N C (L) = ^ clusters of singular energy EJL) = £ e 

SI 

The singular energy density e s i ng (L,W) will then present a complicated finite-size form involving the two ratios 
L/£i and £/£ 

e„„ (1 , l y, fl = 0) = Mpe = |*(|;|) (94, 

where the function $(xi,a;2) should describe the crossover between (i) and (ii) when X2 = 

= 0,x 2 = 0) =1 

$(xi -> +oo,x 2 = 0) (95) 
and should describe the crossover between (ii) and (hi) when x\ = 

$(xi = 0,x 2 ->• +oo) (96) 

x 2 

Note that in usual thermal transitions where there is no droplet exponent 9 = 0, the crossover between (ii) and (hi) 
actually disappears, so that the finite-size behaviors are only governed by the single ratio L/£i as proposed in |58l - [63 |. 
However in the presence of the droplet exponent 9 > for the RFIM, we expect that this simple recipe does not work 
anymore. In particular, the hnite-size behavior exactly at criticality W c (Regime (i) since £i and £ are inhnite) 

e sing (L,W c ) = — d =L 9 - d (97) 

is not connected to the result for W ^ W c in the thermodynamic limit L — > +oo (Regime (hi) since £i and £ are 
finite) 

e smg (L = oo,W ^ W c ) = K = £ e - dc = £— ~ = \W - W c \ 2 (98) 

£i 
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via a single crossover. 

Let us now translate this scaling analysis for the finite-size properties exactly at W c as a function of the external 
field H, to clarify the finite-size properties of the magnetization and of the susceptibility that can be obtained from 
the singularity of the energy density by successive derivation with respect to H. 



B. Finite-size properties at W c as a function of H 




FIG. 8: Finite-size scaling of the magnetization m(H, L) at W c as a function of the external magnetic field H for a = 0.2 
(mean-field region) : (a) L 1//2_CT m(H, L) as a function of the scaling variable HL 1 ^ 2 (Eq. I112[) (b) H ~ 1//3 m(_ff, L) as a function 
of the scaling variable #L 3ct/2 (Eq. [113)1 . 

In section lllBi we have described the mean-field critical exponents for thermodynamic observables at H = 0. Let us 
now describe what happens as a function of the external field H exactly at W c in the thermodynamic limit L — > +oo 

e sing (W c ,H) ~\H\^ +1 
m(W c ,H) 

X(W C ,H) ^\H\^^ (99) 



in terms of the exponent 



Smf = 1 + = 3 (100) 

PMF 



To reproduce the correct divergence of the susceptibility x(W c ,i7) ~ \H\ 2 / 3 via the analog of Eq. [20l the analog of 
Eq. [TO] for the Fourier transform of the connected correlation reads 

C CO nnected{q) ~ |^| 2 /3 + LJ^ ( 101 ) 

so that the exponential decay for H ^ is governed by 



C connected (r) ~ / dge^ -—!— cx e /"> with ^ ~ |JJ|-""i and i/j£> = A (102) 



The the density n(s) of clusters of size s per spin of Eq. l60l becomes 



n ^ = -^+i e ~ SlHli with ™ F = I ( 103 ) 



18 



In a volume l d , the number N t d of clusters has the following singularity in terms of the cut-off s max ~ \H\ 3 of Eq. 

una 

N ld = / d ]Tn( S ) = l d [s max ]-™ F = l d \H\^ M " = l d \H\ 2 (104) 

S 

So the length-scale £1 introduced in Eq. [551 now diverges as 

d[ H) = \H\~i (105) 

The analysis leading to Eq. [MJ is now the same, provided one use i.e. we 

may write after changes of variables to make clearer the dependence upon if 

e sing (L, W c , H>0) = L e H 2 £ fff>/ 2 ^ (106) 

where the function £ (2/1,2/2) satisfies 

£ (yi = 0,2/2-0) =1 

£(j/i -> +00,2/2 = 0) ~ 2/1 

£(j/i = 0,272 ^+00) ~2/^ (107) 

to reproduce the three regimes L <C £1, £1 <C L <C £ and 5«£ 

e sin9 (L,W c ,ff >0) ~L e - d for X « ff- 2/d 
e sin9 (L,W c ,ff >0) ~L e fi 2 for ff ~ 2 / d «L« £ff- 2 / d <- 

e smg (L,W c ,H > 0) ~ H 2-2e/d c = H 4/3 { 0T H- 2 / d "<L (108) 
By differentiation with respect to if, we thus expect the following finite-size behavior for the magnetization 

m(L, W c , if > 0) = L e HM (^^; fff> /2 ^) (109) 

where the function .M (2/1,2/2) satisfies 

X (2/1 =0,2/2 = 0) =1 
M(yi -> +00, 2/2 = 0) ~ 2/1 

■Mfei = 0, 2/2 -> +00) ~ 2/ 2 "^ (110) 
to reproduce the three regimes L <S £1 , £1 L <S £ and £ -C £ 

m(L, W c , if > 0) ~ ^e-d/2 for L ^ jj-2/d 

m(L, W c , if > 0) ~ L e if for ff- 2 / d < L < £H~ 2/d ° 

m(L,W c ,H>0) ~H 1/3 forif" 2/dc <i (111) 

To test the presence of the two crossovers of Eq. Illll for the magnetization, we have plotted our data for the Dyson 
hierarchical model for a = 0.2 in terms of the ratios and L/£ respectively. On Fig. [8] (a), we show the test of 
the scaling form involving 

L 1/2 - a m(L,W c ,H > 0) =Mi (hL^\ (112) 
whereas on Fig. [8] (b), we show the test of the other scaling form involving L/£ 

H- 1/3 m(L,W c ,H > 0) = M 2 (hl¥) (113) 



By another differentiation of the magnetization of Eq. 1 1 1 II with respect to if, one obtains that the crossover at 
£1 actually disappears in the finite-size behavior of the susceptibility (i.e. the regimes (i) and (ii) give the same 
contribution) 

X {L,W C ,H > 0) = L e G(HL d "' 2 ] (114) 
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where the function G(y) satisfies 

G(y = 0) =1 

G(y 2 ->+oo) ~iT^ (115) 

corresponding to the two regimes L <C £ and £ <C i 

VF C , if > 0) ~ L e for L < £H- 2/d " 
X {L 1 W C ,H>Q) ~H- 2/3 lor H- 2/d - <^L (116) 

VII. CONCLUSION 

In this paper, we have discussed the statistics of low-energy excitations and of equilibrium avalanches for the 
zero-temperature finite-disorder critical point of the Random- field Ising model (RFIM). Besides its role is the ther- 
modynamics, the droplet exponent is expected to govern the distribution dl/l 1+e of the linear-size I of low-energy 
excitations or zero-field avalanches. In terms of the number s ~ l df of spins of these excitations of fractal dimension 
df, the power-law distribution thus reads ds/ s 1+e ^ df . In the non-mean-field region d < d c , droplets are compact 
df = d, whereas in the mean- field region d > c? c , droplets have a fractal dimension df = 26 leading to the well-known 
mean-field result ds/s 3 ^ 2 . We have also discussed in details finite-size effects, both in the non- mean- field region 
d < d c where standard finite-size scaling is valid, and in the mean-field region d > d c , where standard finite-size 
scaling breaks down, because the correlation length is not the only relevant length scale. We have explained how to 
adapt the Coniglio's scaling approach to understand the finite size properties of the RFIM for d > d c in terms of the 
two scales £i and £. All expectations have been checked numerically for the Dyson hierarchical version of the RFIM, 
where large systems can be studied with a good statistics via exact recursion, and where the droplet exponent 6 can 
be varied as a function of the parameter a of the effective power-law ferromagnetic coupling. 

Appendix A: Exact recursion for the Dyson hierarchical RFIM 

In this Appendix, we explain how the sequence of ground states that occur as a function of the external field if 
in any given disordered sample can be easily computed by recursion for the Dyson model introduced in Eq. 1281 All 
numerical results presented in the text have been obtained by this method. 

1. Exact recursion for the partition function 

Following the method of [25[ for the ferromagnetic model, we use the Gaussian identity 

e M^) 2 ] = fi r dxe ^ 2 ^{^ *) (A1) 

to rewrite the partition function for a system containing 2^ spins of Hamiltonian given by Eq 1281 

Z 2 N{p;H;{hi,...,h 2 N})= e -m 2 N(H-{h u ...,h 2N }) ( A2 ) 

Si=±l;...,S,jv=±l 



as 



U f + °° 

Z 2 N((3;H;{h l7 ...,h 2 N}) = \j ~ J dxe ~ 



/3x 2 

xZ 2 n-i(/3;H + 2x^J N ;{h 1 ,...,h 2 N-i}) x Z 2 n-i(J3; H + 2x\/J^\ {hus-i+i, H 2 n}) (A3) 
On the right hand-side appear the partition functions of the two half-systems in the modified exterior magnetic field 

H' =H + 2xy/dN. 

The initial condition corresponding to N — and a single spin 2° = 1 reads 

Zi{fi;H;{hi}) = e' 3(H+/ll)Sl = 2 cosh + fci)) (A4) 

Si=±l 
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2. Limit of zero-temperature 



In the limit of zero-temperature where /3 = 1/T — > +00, the partition function of Eq. IA2I becomes dominated by 
the ground-state energy E®n{H\ {hi, h 2 N}) 



Z 2 N{p]H;{hx,...,h 2 N}) ~ e 

,9->+oo 



-/3BS?(if;{/ii,...,h 2jv }) 



(A5) 



In the limit (3 —> +00, the recursion of Eq. IA3I can be thus evaluated via the saddle-point approximation. This 
yields the following recursion for the ground state energy 

E 2N (H; {hi, h 2 N}) 

x 2 + E&-x(H + Ix^lhs; {hi, h 2 N-i}) + E^- X {H + 2x^/1^; {h 2 N-i +x , h 2N }) 



mm/// 



- — — — - — I- E 2 x-i(H'; {hi,..., h 2 N-i}) + E^n-i {H 1 ; {h 2 N-i +1 , h 2 n}) 



(A6) 



The initial condition (Eq. IA4[) reads 



Ei(H;{hi}) = -\H + hi\ 



(A7) 



3. Sequence of ground states as a function of the exterior field H 



It is convenient to characterize each disordered sample of 2 n spins by its 'sequence' of ground states as the ex- 
terior field H is swept from (—00) to (+00) : this sequence contains a certain number (1 + p ma x) of configurations 
{Co, Ci, C Vmasc {, where Co is the configuration where all spins are negative Si — —1 (ground state when H — > —00), 
and where C Pmax is the configuration where all spins are positive S t = 1 (ground state when H — > +00). 

The energy of each configuration C p depends linearly on the exterior field H, with a slope determined by the 
magnetization Mc p = J2i &i 01 ^ ne configuration C p 

E(C p ,H) = -M Cp H + a Cp (A8) 

The value of the field H = Hc p ,c p+1 where the ground state changes from C p to C p +i is simply given by the intersection 
of the corresponding two lines (Eq. IA8|) 

ac p+1 ~ ac P , 

H Cp ,c p+1 = 7 — — r (A9) 

(mc p+1 - run) 



4. Notion of 'no-passing rule' 



The notion of 'no-passing rule' has been first developed for non-equilibrium dynamics concerning charge-density 
waves [63], driven elastic manifolds [64| and the non-equilibrium dynamics of the RFIM [65]. This notion has been 
then extended to the equilibrium of the RFIM [ll|, |66| , where it means that the sequence of ground states that appear 
as a function of the external field H are 'ordered' in magnetization 

M Co = -T < M Cl <M C2 <..< M Cpmai = 2" (A10) 

Since the difference between the magnetizations of two consecutive configurations is bounded from below by the value 
2 that corresponds to a single spin-flip 

- m Cp > 2 (All) 
so that the number p max is bounded by the number of spins 

1 < Pmax < 2" (A12) 

i.e. it always remain very small with respect to the total number 2 2 of possible configurations. 
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5. Recursion on the sequence of ground states as a function of the exterior field H 

Let us now assume that we know the sequences of ground states in H of two independent half-systems of size 2 N ~ 1 , 
and we wish to construct the sequence for the whole system when these two half-systems are coupled. 

We consider the set of pairs of ground states of the two subsystems that exist at a given same exterior field H' : 
let us call I PltP2 the interval of the exterior field H', where the first subsystem has for ground state the configuration 
of energy (Eq. Eg} 

£§L (#') = -m c (u H 1 + a£> (A13) 

(2) 

and where the second subsystem has for ground state the configuration Cp 2 of energy 

£gL (H') = -m c g) H' + a£> (A14) 
The recursion of Eq. IA6I means that in this interval I Pl , P2 of H' , the function that has to be minimized reads 

= ^—77—^ - m rm H' + a W - m (2) R< + a ( 2 ) (A15) 



The minimization over H' 



yields the solution 



Q = d H 4i P1 , P2 {H') = {H ' - H) - m c(1) - m c(B ) (A16) 



ff^(ff) = H + 2J N (m rW +rn r(2) ) (A17) 

that corresponds to the value (Eq IA15|) 



, {Hl(H)) = Jjv(m r (i) + m <2)) 2 - (m ( i> + m c m){H + 2Jjv(m r (i) + m (2 ))) + a£/ + a < - 2) 



Cp 2 



= -Jjv(m „(i) + m„(2>) - (m (i) + m ( 2))H + oi, J + 1 

L'pj ^Pl P2 

= E 2J ,(C = (£7W > CW,H) (A18) 

representing the energy of the global configuration C = (Cpi ,Cp2 ) made of Cp\ and Cp^ for the two sub-systems, 
when the exterior field is H. 



We now need to minimize over all possible intervals I Pl >P2 for H' (Eq. IA6|) 

£ 2 G /(ff) = riun^i, ^ [e 2 «{C = (CW,C«, JJ)] (A19) 

where the minimization is over all the pairs (Cp* , Cp 2 , ) of configurations that are ground-states of the isolated sub- 
systems for some same value H' of the exterior field. 
In practice, we have thus used the following procedure : 

(i) We make the ordered list in the field H' of the pairs (CpP , Cj$ ) of configurations that are ground-states of the 
isolated two sub-systems. Let us we call {Co, C\, Cq max } these 'candidate' configurations of the whole system, and 

compute their energies. For instance if C q = (C p \ , Cp% ), its energy is simply 

E 2 » (C q = (C£> , C$ , H) = -M Cg H + a Cq 

M Cq =M^+M^+a%l+4l 

a Cq = -Jn(M cM + M c<2) ) 2 + ag + c$ (A20) 

Pi P2 PI P2 

(ii) Now these energies are in competition to be the ground state of the whole system at some given exterior field H . 
Since they are ordered in magnetization, one may proceed as follows [IT, . 66] . We start from the two known extremal 
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configurations : the first configuration is Co = (Cq ,Cq ) where all spins are negative (ground-state for H — ¥ — oo) 
and the last configuration is C/ ast = (C^L* , C^. ax ) where all spins are positive (ground-state for iJ — » +00). We 
compute the crossing field H^jast where the two corresponding energies cross. We now compute the energies of all 
intermediate candidates at this crossing field Ho i ast and select the minimal value : the corresponding configuration 
Cq a is then the ground state at Ho^ ast . We may now iterate this procedure : we compute the crossing field Ho Qa 
and find the minimal energy at this field among the candidates < q < q a ; similarly we compute the crossing field 
Hq at iast etc... This method allows to compute the sequence of ground states that really appear as a function of H 
for the whole interval. To make even clearer this procedure, we now describe as an example the first step where two 
systems of one spin are coupled to form a system of 2 spins. 



6. Example with TV = 1 corresponding to 2 1 = 2 spins 

Each subsystem contains only one spin. So the sequence of the ground state of the first subsystem as function 
of the exterior field contains only the two configurations (Cg , C} 1 ^) corresponding to (Si = — I, Si = +1) and the 
energies read (Eq. IA7I and IA8I ) 

E(C { 1] , H') = -m nW H' + a nm = H' + K 

E(C[ 1] ,H') = —m c {\)H' + a a (i) = —H' — hi (A21) 



so that the field H' = H <i) where the ground state changes from to is 



H c w c w = -hi (A22) 
Similarly, the sequence of the second sub-system is described by the parameters 

E(C { 2 \h') = -m n wH' + a„ ( 2) = H' + h 2 

E(C[ 2) ,H') = -m c wH' + a c p = -H> - h 2 
H c (2) c (2) — —I12 (A23) 

Let us assume that hi > h 2 (otherwise exchange the labels 1 and 2). Then H <i) (i> = —hi < —h 2 — H (2) „(2> 

so that the ordered list of candidates for the ground states of the whole system of two spins is (Co = ( );Ci = 

(H — ); C2 = (++)• The energies of these three candidates are 

E(C ,H) = 2H + hi + h 2 - 4 Ji 
E(Ci,H) =-hi + h 2 

E(C 2 ,H) = -2H - hx - h 2 - 4 J x (A24) 
The crossing field Hc ,c 2 where the two energies of the extremal ground states cross reads 

Hc ,C2 = (A25) 

and the corresponding energy reads 

E(C , H Co , C2 ) = E(C 2 ,H Cn . C2 ) = -4Ji (A26) 

We now have to compute the energy of the intermediate candidate C\ at this crossing field 

E(C 1 ,Hc ,c 2 ) = -hi + h 2 (A27) 

to see if it is lower or bigger than the crossing energy of Eq. IA27I : 

(i) if E(Ci, Hc ,c 2 ) < E(Cq,Hc ,c 2 ) = E(C 2 , Hc ,c 2 ): then C\ is indeed the true ground state at the crossing 
field Hc ,c 2 - So the sequence of ground states contains indeed the three configurations (Co, Ci, C 2 ) with the frontiers 
given by the crossing fields Hc ,Ci and Hc 1 ,c 2 - 
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(ii) if E(Ci,Hc ,c 2 ) > E{Co,Hc ,c 2 ) = E{C%, Hc ,c 2 )i then the candidate C\ has to be eliminated. The sequence 

of ground states only contains the two configurations (C = ( ),C 2 = ++) with a frontier given by Hc 0j c 2 - Note 

that the condition for this direct jump between the two ferromagnetic configurations reads 

4Ji > hi -h 2 (> 0) (A28) 

i.e. it is more probable at low disorder as it should. 
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